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The relativistic precession can be quickly inferred from the nonlinear polar orbit equation without 
actually solving it. 
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I. INTRODUCTION 

The precession, and consequently the secular motion 
of the pericenter of a bound orbit, stems from the gen- 
eral relativistic treatment of the motion of a test body in 
the space-time of a spherically symmetric distribution of 
mass which, if also spins, drags space-time around with 
it slightly perturbing the orbit (Lense-Thirring gravito- 
magnetic force), and then by providing a further con- 
tribution to the precession. The study of the former 
phenomenon has been done for various astronomical sce- 
narios. Within the solar system, have been considered 
Earth's satellites?^ the Moor*^^, Mars^ the giant 
planets,±&"^i the Sun and its planets^-22 The 1PN post- 
Newtonian, Schwarzschild-like orbital effects^ recently 
revamped because of several attempts to detect them in 
different natural systems as the galactic environment for 
stellar orbits around the central black holes^r— Still 
in a planetary setting, the exoplanets may constitute a 
fruitful field of study^— 

The geodetic equation in the Schwarzschild space-time 
becomes a Binct-typc differential equation which de- 
scribes in polar form the shape of the orbit of a test 
body under the effect of a force inversely proportional to 
the square of the distance from the origin with added a 
very small inverse fourth power term. There is no diffi- 
culty in principle, using techniques of perturbation the- 
ory, to deduce the precession by solving the equation for 
bound orbits to arbitrary degrees of approximation.— ~— 
Our aim here is to show a new way to obtain the leading 
term of the precession in two moves, consisting of a sub- 
stitution followed by an elementary definite integration 
but, to be able to appreciate its correctness, it needs to 
be opportunely introduced and justified. 



II. THE ANGULAR PERIOD OPERATOR 

We start from the classical equation of the unperturbed 
orbit in polar coordinates^ 

u" + u = A, (1) 

where the primes denote double differentiation with re- 
spect to the angle <f>. The function u = «(</>) is the inverse 
of the radius vector r, while 



where \i is the gravitational parameter of the central body 
(this is strictly correct only when the orbiting body is 
considered to be a test-body) and I is the area constant. 
Equation ([T]) admits u — A as particular integral, and its 
complete solution can be found solving the homogeneous 
equation (u — A)" + (u — A) =0. The solution in the 
perifocal system can be written in the form 



u = A + Ae cos < 



(3) 



which is the equation of the elliptical orbit when e e (0, 1). 
In this setting the angle <f> is the true anomaly while, in 
terms of the elliptical elements A = u(w/2) = l/a(l — e 2 ), 
where a, e are the semimajor axis and the eccentricity. So 
1/A is the semi-latus rectum and determines the size of 
the orbit. The fixed position of pericenter is at <f> = 0, 
where u = 1/ r assumes its maximum value. The solution 
u is periodic of period 2tt, and its average value is 



1 

2^ 



ud(f> = A, 



(4) 



which is precisely the particular integral of Eq. ([T]). It 
corresponds to the circular orbit of a body with the given 
angular momentum, so that we may associate to the el- 
liptic orbit this particular circular orbit. They of course 
share the 27r-periodicity, and this means that, consider- 
ing them together, if we assume that the two radii vectors 
overlap at (j> = 0, the pericenter of the elliptic orbit, they 
will overlap again at the successive pericenter <f> = 2-7T, 
despite from dynamical point of view the respective po- 
lar angles assume over time different values in the other 
positions (apoccntcr excluded) in a way dictated by the 
constancy of the respective areas swept out. But, putting 
aside the time which is extraneous to the mathematics of 
our problem, it is sufficient to think that the two or- 
bits are connected to one another as specified only at 
this special points, and that therefore in passing from 
one pericenter to the successive, both radii vectors ro- 
tate through an angle 2-7T about the origin. This angle 
may be obtained directly from the solution swapping A 
and 2-7T in Eq. so that we may write 



Pu 



ud(f> = 271", 



(5) 



A = 



fJLU 



r- 



>o, 



(2) 



because the definite integration may be considered as the 
application of a linear operator P to the solution u, be- 
ing P a mathematical device acting as a detector which 
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measures the angular period, that is the angle separating 
two successive passages of the body through pcriccnters. 

III. THE RELATIVISTIC PRECESSION 

The general relativistic corrections to orbital dynam- 
ics imply that the classical Eq. ((T|) is substituted by the 
nonlinear polar orbit equation^ 

u" + u r = A + Bui, ( 6 ) 

where B = 3a, being a = 3/i/c 2 the gravitational radius 
of the central attractor, generally a tiny fraction of typi- 
cal orbital dimensions. Without any attempt to try find- 
ing an approximate bound solution, let's handle it in a 
simple, but meaningful way. To do this, wanting to high- 
light a secular effect on the orbit, we put u 2 = u 2 = A 2 
in the last term of Eq. ([6]). This means assuming as 
first approximation to the solution just that circular or- 
bit we associated to the unperturbed elliptic orbit, with 
the shared property with the latter we alluded to before. 
We obtain so a linear equation containing a constant per- 
turbation 

u" + u r = A r , A r = A + A 2 B, (7) 

whose solution in the perifocal system may be written in 
the form 

u r = A r + A r e cos 4> — u + A 2 B + A 2 B e cos 4> 
= u + Su, (8) 

namely as the sum of the unperturbed elliptic solution 
with added a perturbation term factored by 5. It is 
worthy noting that e is unaffected by the perturbation, so 
it is still completely arbitrary in its interval of definition. 
As we consider Eq. J7| as a perturbed version of Eq. (TTJ) , 
it is legitimate to apply to the solution ((8|) the operator 
P, the detector of the angular period, and thus we get 

Pu r =Pu + PSu = 2ir + 2irAB, (9) 

which also represents the lowest-order contribution to the 
precession of the circular orbit per revolution that occurs 
in the direction of increasing true anomaly. In substance, 
the relativistic term makes a circumference of the circular 
orbit shorter by AB which in turn makes a deficit of an- 
gle of rotation of 2ttAB. So, to complete one revolution 
along the orbit of radius 1/A for which is calibrated the 
operator P, the orbiting body needs to rotate this same 
extra angle and this in turn is reflected in the pericenter 
shift of the associate elliptic orbit because of the compul- 
sory coincidence of the angular positions of the two radii 
vectors at this special point. In conclusion, expressing 
A in terms of the elliptical elements of the orbit, from 



Eq. ([9]) it follows the Einstein's precession formula 
6na 

Pu r = 2ir+— 10 

a(l — e z J 

The essential feature of the approach followed (a par- 
ticular averaging method) is that, without the need of 
using any specific perturbation equation, it employs the 
circular motion as a test motion for which, thorough the 
operator P, is easy to obtain the leading precession term, 
which is pari passu transferred to the elliptical orbit. 
When applied to other types of inverse-power perturbing 
central forces (not necessarily of academic interest) this 
method gives the dimcnsionally correct algebraic form of 
the precessional term but off target by a small numerical 
factor, so underestimating the proper orbital precession 
determined with other methods4£~— We deduce that in 
these instances the procedure cannot fully capture the 
real variation of A on the sole basis of the unperturbed 
averaged value of u employed in the approximation, and 
this is perhaps due to the fact that the constant A ap- 
pears linearly after the action of the operator P only with 
the occurrence of a quadratic nonlinearity as in Eq. (|6|) , 
and we were lucky enough that in the most important 
case of physical interest it works well. In the other cases, 
the method should be refined. 



IV. CONCLUSIONS 

It would be interesting to see whether it can similarly 
treat another case of possible relativistic precession, the 
Lense Thirring drag^i to which we dedicate here only 
a fleeting mention. This effect closely parallels the way 
in which a rotating electrically charged body generates 
magnetism. The analogy has helped to organize the un- 
derstanding of the phenomenon and to determine pre- 
dictions experimentally testable. The Gravity Probe B 
satellit o 58 ' 59 achieved an accuracy within 19 per cent of 
the expected orbit change; other satellites got a rather 
similar level of accuracy. Researchers hope to achieve 
1 per cent with LARES.— ~— In this phenomenological 
approach the forces in play are non-central, so that our 
method is not directly applicable, but it would not be 
difficult to conceive a model in an appropriate plane of 
symmetry which allows to simply deduce, with suitable 
geometrical and physical insight to justify the mathemat- 
ics employed, the leading precessional term with a mini- 
mum effort. We suggest also to use the method in some 
non-Newtonian model of gravity, as the Yukawa-like fifth 
forced Here we have no problems of non-centrality, so it 
is easy write down the relative polar equation and hence 
the precession term, but after that assumptions must be 
made on the figures to be assigned to the free parame- 
ters in such a way as to obtain testable results for realistic 
models of astrophysical systems. 
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